A fundamental limit for integrated atom optics with Bose-Einstein condensates 
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The dynamical response of an atomic Bose-Einstein condensate manipulated by an integrated 
atom optics device such as a microtrap or a microfabricated waveguide is studied. We show that 
when the miniaturization of the device enforces a sufficiently high condensate density, three-body 
interactions lead to a spatial modulational instability that results in a fundamental limit on the 
coherent manipulation of Bose-Einstein condensates. 
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I. INTRODUCTION 

The central idea of integrated atom optics is to minia- 
turize the size of active and passive atom optical com- 
ponents such as atom lasers, wave guides, beam split- 
ters, and interferometers, and to combine them in in- 
tegrated devices [llllllllIllllSli|lS|llE2. 
Potential future appHcations include hand-carried high- 
precision measurement devices such as gravimeters and 
gyroscopes. Integrated atom-optical chips may also prove 
useful to control neutral atoms in integrated and scal- 
able microtrap arrays for quantum information process- 
ing. For example, microfabricated atom optical compo- 
nents can be fabricated with current-carrying wires or 
microstructured surfaces, combined with homogeneous 
magnetic bias fields and periodically magnetized sub- 
strates. Optical fields or electric fields may also be in- 
volved in the hybrid construction of such devices. 

However, the future of integrated atom optics still de- 
pends on finding solutions of a number of both techni- 
cal and fundamental issues. For example, it will be es- 
sential to maintain the coherence of the matter waves 
to be manipulated, a step still to be demonstrated in 
waveguide-based beam splitters. Atoms confined in wire- 
based microtraps on atom chips face a relatively noisy 
environment, due e.g. to the heating from the chip sub- 
strate, which can cause the decoherence of atomic mat- 
ter wave III addition, fluctuating currents in micro- 
wires may perturb the atoms when they are brought near 
its surface: Recent experiments have revealed the frag- 
mentation of atomic Bose-Einstein condensates (BECs) 
confined in wire-based microtraps when the atom-surface 
distance is reduced to micrometer scale jl^ JJ^ , and 
the experimental data seems to support the argument 
that the fiuctuating currents in the wire and corruga- 
tions due to imperfect microfabrication are responsible 
for that fragmentation. 

While such experiments illustrate the technical limita- 
tions for the coherent control of atoms in the microtraps, 
it can be hoped that these will be eliminated by bet- 
ter engineering solutions. In addition, though, it is also 
important to understand the fundamental limitations to 
integrated atom optics. For example, the strong com- 
pression typical of microtraps or waveguides can signifi- 



cantly increase the density of the atomic BECs, thereby 
enhancing inelastic collisions and resulting in a funda- 
mental source of decoherence that reduces the lifetime of 
the confined BECs (lOj . 

In this paper, we discuss an additional effect that lim- 
its the coherent control of atomic BECs in microtraps or 
microfabricated waveguides. Specifically, we show that 
when the miniaturization of the device results in a suffi- 
ciently high condensate density, three-body interactions 
lead to a spatial modulational instability that leads to 
the fragmentation of the condensate, thereby imposing 
a fundamental limit on the densities that can be stably 
propagated and manipulated in atomic waveguides. 

The paper is organized as follows: Section II introduces 
our model and establishes the notation. Section HI uses a 
linear stability analysis to show analytically the existence 
of a modulational instability in the condensate dynamics 
when three-body collisions become important. Section 
IV presents the results of selected numerical simulations 
that confirm the analytical predictions and show the on- 
set of condensate fragmentation. Section V expands our 
analysis and presents arguments that reinforce the con- 
clusions of the numerical simulations and argue that the 
modulational instability eventually leads to a collapse of 
the condensate. Finally, section VI is a summary and 
outlook. 



II. THE MODEL 

For concreteness, we consider in this paper an atomic 
waveguide resulting from the combination of a current- 
carrying wire (for instance) electroplated on a substrate, 
and a static, homogeneous bias magnetic field Sbias 01 
perpendicular to the conductor. Assuming a linear wire, 
the distance d of the waveguide to the chip surface is 



\2ttJ Sbias ' 



(1) 



where fio is the vacuum permeability and / the current 
through the wire. Since the guide height can be reduced 
using smaller currents or stronger bias field, the setup is 
ideal for the requirement of miniaturization in integrated 
atom optics. 
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The radial oscillation frequency of this trap is 
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where /^s is the Bohr magneton, gp and mp the gyro- 
magnetic ratio and magnetic moment of the hyperfine 
state of the alkali atoms of mass m being trapped, and 
we have assumed that the axial confinement of the atoms 
is achieved by a loffe-Pritchard-type trap that provides 
a field i?o in the center of the trap. 

Since the radial oscillation frequency ujr scales as 
^bias/^' ^"^y miniaturization of the trap achieved through 
a reduction of d results in an increase of Wr, and hence 
in the density of the confined condensate. The closer to 
the chip surface the condensate is brought, the higher 
its density. In some of the current experiments, the con- 
densate densities can already reach values in excess of 
1 X 10^^ cm~'^. For such high densities, the effects of 
three-body interactions on the dynamics of BECs are ex- 
pected to become important. Our goal in this paper is to 
understand their impact on the dynamics of the guided 
condensate. 

Our starting point is the three-dimensional Gross- 
Pitaevskii equation for the wave function 5'(r,t) of a 
trapped condensate, generalized to include three-body 
interactions. 



(3) 



Here fi is the chemical potential, (72 and ^3 measure the 
strength of the two-body and three-body interactions, 
respectively, \l/(r, t) is normalized to unity and N is the 
number of atoms in the condensate. We have also intro- 
duced the dimensionless ratio 



(4) 



which determines the strength of the transverse confine- 
ment relative to the longitudinal trapping. (We allow 
£_{t) to be time-dependent to include possible time varia- 
tions of bias magnetic field Bbias or of the current /, and 
hence of w^.) 

In the following, we assume that the microtrap is op- 
erating in the regime of tight transverse confinement, 
^{t) » 1 and also that the transverse trapping poten- 
tial is much larger than the interatomic interaction. We 
further assume that the time variation of LOrit) is slow 
enough that the transverse profile of the BEG adiabati- 
cally follows the ground transverse oscillator state 



(5) 



We can then proceed by approximating the condensate 
wave function ^'(r, i) as 



(6) 



where it is understood that the radial wave function 
(j)o{r,t) contains only slow time variations, and the en- 
velope function u{z,t) incorporates all other time varia- 
tions. We then obtain the approximate equation for the 
envelope function 



^ du 
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-%2,2(t)|u|'M + %3,.Wk|*U, (7) 



where the reduced interatomic interaction coefficients 
g2,z and 53, ^ are 

92At) = 92N [ 2nrdr\Mr,tr=(^]m: 



93,z{t) = gaN^ / 2TTrdr\Mr,t)\ 



93N' 
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with I2 = ^Jh/niLOz- From these definitions we see that 
the reduced two-body interaction coefficient varies lin- 
early with the transverse oscillator frequency or the ra- 
tio ^(t) = ujt{t)/u}z, whereas the three-body coefficient 
scales as (i) . This implies that the relative significance 
of two- and three-body interactions depends on the level 
of transverse confinement. 

For alkali atoms, the two-body interaction coefficient 
has fully been determined by experiments. Although in- 
formation for three-body interactions is still sketchy, re- 
cent theoretical studies have provided a way to estimate 
the three-body interaction coefficient gs 0, \l5, 16, ij, 
IT^ . For rubidium atoms, it may be expected to be nega- 
tive, with 1^3 1 of the order of IQ-^** - 10"" cm^/s. With 
this in mind, we focus our discussion on atoms with re- 
pulsive two-body interactions (32 > 0) and attractive 
three-body interactions (f/3 < 0). 



III. MODULATIONAL INSTABILITY 

To expose the basic modulational instability that arises 
at high densities we first consider the case where the lon- 
gitudinal trapping potential is switched off, the trans- 
verse trapping is fixed, so that ^ >> 1 is constant, and 
the condensate has a homogeneous longitudinal distribu- 
tion uq over a length L. The chemical potential is then 
M = ^92,z\uo\'^ - fl\93,z\\uo\'^- 

We are interested in the condensate dynamics in the 
presence of the attractive three-body interaction. For 
this purpose, we first carry out a linear stability analysis 
by setting 



with lr{t) — \/h/mujr{t). 



u{z,t) = [uQ + Su{z,t)]e-"''^'', 



(9) 
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and looking at longitudinal excitations of the form 

k 

Substituting Eqs. and ifTUI) into Eq. we obtain 

VkUk 

-VkUk 



h 52 



where we have defined 



(10) 
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Considering plane- wave excitations Uk = i^/L) exp(iA;z) 
and Vk — (1/i) exp(— Eq. ifTUIl gives the excitation 
frequencies 



i^k 



2to V 2to 
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which reduce to the familiar Bogoliubov spectrum in the 
absence of three-body excitations. In that limit, the exci- 
tation energies are always positive and the condensate is 
stable. However, the situation changes completely when 
92, z < 2|(73^z||moP, or 77 < 0, in which case the excitation 
energies can become imaginary. Physically, such complex 
frequencies imply that spatially modulated perturbations 
can grow with time, that is, a modulational instability 
may occur. The maximum gain occurs for the wave num- 
ber k such that Tik'^/2m = |?7|, which yields the spatial pe- 
riod of the most unstable wave as L„i = ^J^n'^hj {m\T\\). 
Provided that the length of the BEC satisfies L > L,,,,, 
it will then become modulationally unstable due to the 
attractive three-body interactions. As a result, a spa- 
tial inhomogeneity will develop in the condensate density 
profile along the z-axis. Three-body interactions impose 
a fundamental limit on the linear densities that can be 
stably propagated and manipulated in the waveguide. 

Consider a *^Rb BEC as an example. The two- 
body interaction coefficient of ^^Rb is (?2 — 'iiT'has/m ^ 
4.95 X 10""'^^ cm'^/s. The theoretical estimate of the three- 
body interaction coefficient is Igs] ^ 10"^^ — 10"^^ cm^/s. 
The onset of the modulational instability occurs then at 
IssIpo ~ 52/2, with po being the three-dimensional den- 
sity of the condensate. This yields a threshold density 
of the order of 10^^ — 10^^ cm~^ beyond which the con- 
densate becomes dynamically unstable. We remark that 
such densities might have already been achieved in recent 
experiments. The corresponding modulational length is 
on the order of a few microns. 



IV. NUMERICAL SIMULATIONS 

In this Section we present numerical simulations that 
illustrate the effect of the modulational instability on 



the condensate dynamics. In particular, we consider the 
case the magnitude of the bias magnetic field -Bbias is de- 
creased to lower a magnetically trapped gas towards the 
atom chip surface. We assume that the initial trapped 
gas is stable against modulational instabilities. As the 
condensate approaches the surface of the chip, its lin- 
ear density increases concomitantly with the increasing 
transverse trapping frequency. The gas can then become 
modulationally unstable. Our numerics give illustrative 
examples of the development of this instability. 

For convenience of numerical calculation, we intro- 
duce the dimcnsionless variables t — ujzt,C, ~ z/l^, and 
*i'(C, ''') = a/U"", to yield 



2 1 (j) 



1 

M{rm^^ + M\r)\^^^, (13) 



where 

P2 = 



<0, /33 = 



Co- (14) 



Here we have written 



^(t) = Wr(r)/u;3 = Co/(t"), 

^0 = '-^r{0)/ujz being the initial ratio of the transverse and 
longitudinal trapping frequencies. The function /(r) in- 
corporates the time variation of the magnetic bias field. 
We choose ^0 such that three-body interactions are ini- 
tially negligible for the trapped gas and the condensate 
is initially stable. As it is lowered the transverse trap- 
ping frequency and linear density increase, and so /(t) 
increases from unity. We note from Eq. (|13|l that the 
two-body interaction term is proportional to /(r) and 
the three-body term is proportional to /^(t), so that as 
the BEC approaches the surface the importance of three- 
body interactions increases, as we have seen. 

We take parameters appropriate to a ^^Rb BEC com- 
posed of = 10"' atoms, m = 1.44 x 10^^^ kg, 32 = 
4.95 X 10"" cm^/s, = -4 x lO^^^ cmVs, with a 
longitudinal trap frequency uj^ = 27r x 14 rad/s, giv- 
ing Iz = 2.8 /Ltm [23. With these parameters we find 
02 = 375 and j/Jaj = 270. We choose the initial trans- 
verse trapping frequency such that ^0 = ^riO)/^z = 10, 
and assume a linear ramp of ^^(t), /(t) = 1 + ar, with 
a = 1.5. The initial macroscopic wave function $(^,0) 
is chosen as the ground state xiC) of the trap, neglect- 
ing three-body interactions. An approximation to the 
ground state is obtained in the Thomas-Fermi approxi- 
mation by setting ^(C, t) = exp(— i7r)x(C) and imposing 
the normalization condition /rfClx(C)P — 1 to yield 



lx(C)l' 



1 

1 m 
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(15) 
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FIG. 1: Scaled density |*&(Ci''")l^ versus and r for the case 
with no ramping of the transverse trapping frequency {a — 0). 

with 9{y) the Heaviside step function. Our condition that 
the transverse trapping frequency varies slowly compared 
to the longitudinal dynamics requires 7 >> a, which is 
satisfied by our chosen parameters. 

Figure n shows the evolution of the scaled linear den- 
sity |$(C, t)P versus ( and r in the absence of ramping of 
the transverse trapping frequency (a = 0). The density 
profile remains intact in time, except for a small mod- 
ulation that arises from the fact that the initial macro- 
scopic wave function is the ground state without three- 
body interactions, whereas the numerical simulation in 
Fig. ^ incorporates three-body interactions. This re- 
sults in some small rearrangement of the density profile. 
The robustness of the density profile to the inclusion of 
the three-body interactions shows that the initial state is 
stable and the three-body effects are not important for 
these initial parameters. This also justifies our use of this 
particular initial wave function. 

This result should be compared and contrasted with 
Fig. 12 which shows the evolution of the scaled lin- 
ear density |$(C,r)p versus r for r = ^ 7.5 for 
a situation where the transverse trapping frequency is 
ramped (a — 1.5). As the transverse trapping frequency 
5(r) = uj{t)/llJz = 10 X (1 -I- 1.5t) increases in time, the 
initial effect is a broadening of condensate width. This 
is a consequence of the fact that the dominant repul- 
sive two-body interactions are becoming stronger. The 
density undulations seen in Fig. |21 arise from the inter- 
play between the attractive harmonic trapping and the 
increasing two-body repulsion. 

As time increases further, the role of the attractive 
three-body interactions starts to become important, with 
the onset of the modulational instability kicking for r > 
7.5. Figure |31 shows the development of the instability 
around the center of the trap at the dimensionlcss times 
T =7.64 (dotted line), r =7.68 (dashed line), and t=7.72 
(solid line) . A spatial modulation of the condensate den- 
sity of period L/l^ — 0.8 is clearly seen to be developing, 
resulting in a fragmentation of the density profile of the 
trapped gas. 

One can compare this period with the analytical pre- 
dictions of Sec. Ill by treating the center of the BEG as 
approximately homogeneous. Converting to dimension- 
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FIG. 2: Scaled density |"l>((;^,r)p versus ( and r for the case 
with ramping of the transverse trapping frequency (a — 1.5). 
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FIG. 3: Scaled density \^{C,r)f versus ( for times r =7.64 
(dotted line), r=7.68 (dashed line), and r=7.72 (solid line). 

less units, the period of the most unstable wave vector 
becomes 



h ^ ^ mm..n,-2\P3\fl,,nl\' ^ ' 

where /max = 13 is the value of /(t = 8) for times in the 
vicinity of where the instability develops, and ric ~ 0.06 
is the value of the scaled linear density at the trap center. 
From this we obtain Lm/lz ~ 0.74, in reasonable agree- 
ment with the period 0.8 between the peaks in Fig. O 
We remark that the threshold for the modulational insta- 
bility occurs when the denominator in Eq. H16|l vanishes, 

that is, 1^2/max'T-c " 2 1/33 | /max"-c I = 0- 

For times r > 8 the main density peaks generated by 
the modulational instability in Fig. |31 continue grow- 
ing, and the numerical simulations quickly break down 
thereafter. Indeed, the numerical solutions indicate that 
the density peaks are undergoing a collapse whereby the 
density locally increases towards infinity within a finite 
time. AbduUaev et. al. |2lj | have previously shown 
that collapse is possible in a trapped gas with repul- 
sive two-body interactions. This collapse is analogous to 
the corresponding case for two- or three-dimensional con- 
densates, where a collapse occurs for condensates large 
enough that the attractive two-body interactions over- 
come the "quantum pressure" associated with the trap 
ground state. By contrast, the present system is effec- 
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tively one-dimensional, (in but with a combination 
of repulsive two-body interactions and attractive three- 
body interactions, and it is the competition between 
these two mechanisms that leads to condensate collapse 
for sufficiently high densities. This places a fundamen- 
tal limitation on the linear densities that can be stably 
sustained in atom microtraps and waveguides. 

In practice, the predicted collapse will be regularized 
by some additional physical process, such as three-body 
losses as in the case of higher-dimensional collapse in the 
Bose-Novae. We conjecture that a form of soliton turbu- 
lence will ensue |22Ll23|| . In addition, as the linear density 
continues to grow the mean-field energy will eventually 
exceed the transverse mode energy, in which case mul- 
tiple transverse modes will be excited thereby arresting 
the collapse. 

V. COLLAPSE AND MODULATIONAL 
INSTABILITY 

In this section, we present a discussion that reinforces 
the argument that the modulational instability eventu- 



ally leads to condensate collapse. We proceed by reex- 
amining the generalized Gross-Pitaevskii Eq. H13|) , which 
can be interpreted as resulting from the effective poten- 
tial 

Consider a normalized Gaussian approximate trial solu- 
tion of width I 

e-i'/^'^'\ (18) 

which relates the peak density |<&(0)p at trap center to I 
via 

ne = |<i>(0)p^l/V^. (19) 

Near the center, and evaluating the potential for param- 
eters /max: in the vicinity of where the modulational 
instability develops, we have 



UeffiC) ~ (/32/maxnc " l/^S |/n' ax"^) + " /^27r/,„ax"c + 2 j/Jg 1 7r/,Lx^c ) +- (20) 



We are interested in the under-braced term, which con- 
trols the dominant confinement properties due to the 
combined effects of the linear trap and two- and three- 
body interactions. At densities low enough for two-body 
interactions to dominate over three-body processes, the 
two-body collisions oppose the focusing of the trap, and 
their balance produces a stable and confined solution. In 
contrast, when attractive three-body interactions domi- 
nate, both the trap and three-body interactions work to 
compress the density, causing a strong focusing of the 
density distribution and collapse. We therefore antici- 
pate that the cross-over from stability to collapse will 
occur when the two- and three-body focusing terms in 
the under-braced term cancel each other, 

/32^/maxnc " 2|/33 k/^^^ax^c = 0. (21) 

This coincides precisely with the modulational instability 
threshold from Eq. (|16|l . Thus, the appearance of a mod- 



ulational instability also signals that the condensate will 
collapse if the trapped gas is too close to the atom chip 
surface. This conclusion is also in agreement with the 
results of AbduUaev et. al. regarding the stability 
of trapped gases. We remark that although we have con- 
sidered a magnetically trapped gas, the same conclusions 
clearly apply to an optically trapped gas at sufficiently 
high densities. 

We can also use a variational approach to treat the col- 
lapse of the condensate more rigorously. Take Eq. H18() 
as the variational wave function with the width I being 
the variational parameter. The energy functional associ- 
ated with the dimensionless Schrodinger Eq. (|13|) is then 
given by 



E{1) 



4 



dz $(C) 

A 
I 



15^ 



B 



(22) 
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where we used Eq. A = /?2//(2\/2^) and B = 

mp/i^v^TT) - 1/4. 

From Eq. we have that E{1) —oo as I 

provided that B > 0. This means that under such a 
condition, the condensate is not stable against coUapse. 
However the system can still become metastable if E{1) 
possesses a local minimum at finite I. This is analogous 
to the metastability of an attractive condensate with a 
sufficiently small number of atoms 19]. The condition 
for the existence of a local minimum is dE{l)/dl — and 
(pE{l)/(Pl > 0, which yields 

F{1) = - Al + 2B = 0, 



It can be shown that the above conditions are satisfied if 
and only if 

^f^)=-f^V--^<0- (23) 
\3A ) 2 yiAJ 3 ^ ' 

The local minimum vanishes when the above inequality 
becomes an equality. For the parameters of our numerical 
calculations this happens at t « 8.4 which is in reason- 
able agreement with the onset of collapse t > 8 found 
in the simulations. This not only further establishes the 
connection between the modulational instability and con- 
densate collapse, but also provides the condition for the 
threshold of the instability through Eq. 123f) . 



two-body collisions and attractive three-body collisions, 
is reminiscent of the physics underlying the collapse of 
condensates with attractive two-body interactions, ex- 
cept that in that latter case, the competition is between 
the effect of collisions ant the quantum pressure associ- 
ated with the trap ground state. In both cases, though, 
the competition leads to a fundamental limit on the size 
of condensates that can be stably trapped and coherently 
manipulated, since the modulational instability and col- 
lapse correspond to a form of spatial fragmentation which 
will persist for trapped gases even though current tech- 
nical sources of noise are eliminated. 

In addition to predicting a fundamental limit in in- 
tegrated atom optics, the present study also opens up 
intriguing new directions of investigation. For instance, 
we recall that the use of Feshbach resonances to switch 
the sign of the scattering length from positive to nega- 
tive has lead to the discovery of fascinating dynamical 
effects such as Bose-Novae in three-dimensional conden- 
sates. Similar studies, but in one dimension, should now 
become possible simply by modulating either the bias 
field -Bbias or the current / in wire microtraps. In prac- 
tice, the predicted collapse, and the subsequent dynamics 
of the collapsing filaments, will be regularized by some 
additional physical process, such as three-body losses as 
in the case of higher-dimensional collapse in the Bose- 
Novae, or excitation of higher-order transverse modes as 
the mean- field energy rises due to collapse. In future 
work we shall investigate the detailed dynamics of this 
one-dimensional Bose-Novae beyond the initial growth 
of the modulational instability considered here. 



VI. SUMMARY AND CONCLUSIONS 



Generic microtraps and atomic waveguides based on 
chip technology are characterized by an increased level 
of transverse confinement, and hence increased atomic 
density, as their distance from the chip is reduced. As a 
result, three-body collisions become increasingly impor- 
tant. We have shown that for attractive collisions they 
can lead to a modulational instability in the dynamics of 
the condensate, and its eventual collapse, in agreement 
with the prior work of AbduUaev et. al. 21] on the sta- 
bility of trapped gases. The physics underlying this col- 
lapse, which involves the competition between repulsive 
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